Large-scale surface reconstruction energetics of Pt(100) and Au(100) by all-electron 
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The low-index surfaces of Au and Pt all tend to reconstruct, a fact that is of key importance 
in many nanostructure, catalytic, and electrochemical applications. Remarkably, some significant 
questions regarding their structural energies remain even today, in particular for the large-scale 
quasihexagonal reconstructed (100) surfaces: Rather dissimilar reconstruction energies for Au and 
Pt in available experiments, and experiment and theory do not match for Pt. We here show by all- 
electron density-functional theory that only large enough "(5xN)" approximant supercells capture 
the qualitative reconstruction energy trend between Au(100) and Pt(100), in contrast to what is 
often done in the theoretical literature. Their magnitudes are then in fact similar, and closer to 
the measured value for Pt(100); our calculations achieve excellent agreement with known geometric 
characteristics and provide direct evidence for the electronic reconstruction driving force. 

PACS numbers: 68.35.B-,68.35.Md,71.15.Mb,73.20.At 



C/5 



C 

O 

i o i 

(N 
> 

oo 

o 
o 



X 



The late 5rf transition metals Pt and Au and their 
surfaces are of central importance for a wide variety of 
applications, including catalysis, electrochemistry, sub- 
strates for nanostructure creation and characterization, 
etc. (some recent examples are discussed as Refs. [H-|2lj|. 
below). Their low-index surfaces (100), (110) and (111) 
are paradigm systems to understand catalytic and elec- 
trochemical processes at the atomic scale, to the point 
that only textbooks and handbooks (e.g., Refs. P, 0) 
attempt somewhat comprehensive overviews. Au(lll), 
(100) and Pt(100) are ubiquitous as substrates for nanos- 
tructure creation and manipulation, both due to their 
relative inertness and due to the fact that they tend 
to form well-defined, large-scale quasihexagonal surface 

but also in 



2J-|27|— not just 



m vacuo 



reconstructions l 

solutions, under electrochemical adsorption^ and epi- 
taxial growth conditions. 0, Q The reconstructions im- 
ply a significant rearrangement of surface atoms, and a 
creation of qualitatively distinct surface areas [2^, [29[ 
[fee / hep stacking for (111), steeper vs. flatter ridges 
for (100)]. These local structure variations of Au and 
Pt (100) and (111) surfaces are the key ingredient for 
subsequent nanostructure growth and processes on to 
e.g ., by way of a preferential nucleation of deposits ( 
Il0| | and refs. therein). To give some recent examples 
where this atomic structure is important: the forma- 
tion of magnetic structures [ill], self- assembled nanos- 
tructured arrays of (bio)molecules,[l2| magnetic dots[H 
or supr amolecules . (1 3l 1 141 ] the structural and electronic 
underpinnings of catalytic reactivity, 0, [|| [TBI on-surface 
synthesis of well-separated molecular wires [16[ for trans- 
port studies, 16, [ljj well-defined patterns of molecular 
switches. [Til 19 1 or atomic-scale visualization of kinet- 
ics and dynamics of surface processes. 0, [H, [H 21 1 The 
balance of reconstruction/deconstruction can be essen- 
tial for chemical processes ([!, HH and refs. therein); for 



instance, the peak of catalytic activity of Pd/Au(100) oc- 
curs at a coverage of 0.07 ML, [30] precisely in the range 
where one would expect reconstruction/deconstruction 



to occur. 
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In order to obtain a truly predictive theoretical under- 
standing of these phenomena, one would ideally like to 
use current first-principles theory for "realistic enough" 
unit cell sizes. However, this endeavour faces a significant 
obstacle: The sheer surface reconstruction size necessi- 
tates that even very recent studies employ instead much 
smaller model supercells. [l^, 13 1 This obstacle is most 



pronounced for the (100) surfaces, whose quasihexago- 
nal superstructures are almost uniformly contracted com- 
pared to bulk (111) planes [25% more atoms per area 
in the "hex" layer than in a bulk-like (100) plane] and 
are, strictly, incommensurate with the underlying square 
substrate. [28, [29L I32H36T ] In vacuo, they are additionally 
slightlyrotated (0.5-1°) against the nominal (100) surface 
rows, 32J, 13514381 ] but the exact "hex" layer geometry and 



surface conditions are interrelated: On ste ppe d surfaces, 
periodicities change, [39j and at high T[35j, |36j, |38| or un- 
der catalytic conditions, @ unrotated "hex" planes are 
observed. Again, this would not be a grave obstacle for 
either current experiments or theory if suitably accurate, 
small periodic approximants for the full surface structure 
could be found. However, it turns out that the smallest 
approximants which cover all key features (2-dimensional 
lateral contraction, rotation, differently buckled qualita- 
tive surface areas) are already significant in size: Roughly 
speaking, "(5x7V)",[lj [U Si S3 where 20-30, i.e., 
>100 atoms to be considered in each layer. 

The primary point of the present paper is to demon- 
strate, from first principles, that indeed large-scale 
"(5xiV)", and not smaller approximants to the Au(100) 
and Pt(100) quasihexagonal reconstructions are essen- 
tial to obtain the qualitatively correct surface structural 



2 



energetics from first principles. Our results are based 
on a numerically converged all-electron description, us- 
ing density-functional theory (DFT) in the local-density 
(LDA) and PBE generalized gradient approximation. 
While computationally challenging, such a description is 
now feasible. In addition, it provides a key benchmark for 
any computationally cheaper, more approximate meth- 
ods to be used. 

Concerning the reconstruction energy and driving 
force, some significant qualitative lessons can already 
be learned from simple (lxl) surface models. Using 
DFT-LDA, Takeuchi, Chan, and Ho (TCH) [H| demon- 
strated that free-standing Au(lll) planes achieve a large 
energy lowering when uniformly contracted to the in- 
teratomic spacing of the Au(100) "hex" reconstruction. 
Using a Frenkel-Kontorova model, TCH found a net 
surface energy gain through reconstruction, albeit very 
small: Ai^hox <0.01 eV per unit area of the (100) sur- 
face [eV/lxl]. This is not far from an experimental, 
electrochemical estimate (0.02 eV/lxl j42[), but, puz- 
zlingly, smaller by a factor «6 than a careful assess- 
ment of the reconstruction energy of Pt(100) in vacuo, 
0.12 eV/lxl [H] (an earlier analysis placed Pt(100) even 
higher, at 0.21 eV/lxl HSU). Adding to the con- 
fusion is the fact that theoretical estimates for Pt(100) 
based on small (5x1) approximant cells in vacuo yield 



only half the experimental value (0.05-0.07 eV/lxl 46- 
48|). For Au(100), (5x l)-based theoretical estimates|49- 
51] are approximately in line with the electrochemical 
experiments. Since Au(100) and Pt(lOO) do behave sim- 
ilarly also regarding thermal stability, 35, 36] one might 
consider the experimental reconstruction energy estimate 
for Pt(100) an outlier — if it were not for the fact that it is 
this experiment which pertains to the full reconstruction 
in the vacuum environment, and not to an electrochem- 
ical environment, or to a restricted unit cell. Remark- 
ably, our calculations of the full (5xN) reconstruction 
energy resolve this puzzle in favor of the Pt(100) exper- 
iment, and show that both the experimental and earlier 



51| underestimate the 



theoretical values for Au(100)|4E 
reconstruction energy. 

Figure [1] shows the unit cells of the specific approxi- 
mants that are important for this work. In each case, 
the hexagonal top layer (circles) is slightly distorted to 
form a commensurate coincidence lattice (arrows) with 
the underlying square substrate (crosses). Sub-figure (a) 
shows the popular (5x1) approximant. Compared to 
an ideal hexagonal layer, the top layer is compressed by 
4 % in one dimension but not in the other. To account 
for the reconstruction in both dimensions, larger (5xN)- 
type approximants (N w 20-30) are needed. In subfig- 
ure (b), the hexagonal plane is thus compressed in both 
dimensions, and somewhat distorted so that its close- 
packed rows are additionally rotated by ~ 1° against 
the substrate. Finally, subfigure (c) shows a closely re- 
lated (5xJV)-type approximant, with a rotation of « 1° 
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FIG. 1: Structure and distortion of the quasihexagonal layer. 
Circles: topmost ("hex") layer, crosses: second (square) layer. 
Unit vectors of the coincidence lattice of both planes are 
shown in red. (a) (5x1) and (b) (5x20) superstructure. The 
"hex" layer is compressed by 4% wrt. bulk (111) in both di- 
mensions, and rotated by 1° in one dimension, (c) (5x20)-like 
superstructure, rotated in both dimensions. 



in both dimensions. We model the "hex" layers as part 
of five-layer slabs (two layers relaxed, three kept fixed in 
bulk positions) that are separated by more than 25 A of 
vacuum. The largest slabs (N=40) thus comprise 1046 
atoms, 446 of which are fully relaxed. We emphasize that 
this number refers to heavy elements, no pseudoization 
is employed, metallic systems are not amenable to cur- 
rent O(N) methods, and that especially transition metal 
slabs can be subject to serious self-consistency instabil- 
ities (charge sloshing). 52] On the methodological side, 
our work therefore reflects significant efforts towards scal- 
ability and efficiency (including a real-space version of 
Kerker preconditioning 531) in our accurate "DFT and 
beyond" code FHI-aims. [531 H3| which was used through- 
out this work on massively parallel hardware (IBM's 
BlueGene/P) and on a recent, infiniband-connected Sun 
Microsystems Linux cluster. 

Compared to a (lxl) layer, a (5xN) plane (N >1) 
contains (iV+5+1) additional atoms in the unit cell. In 
terms of total energies for individual surface slabs, £J slab ; 
and the total energy per atom in the bulk, £ ,a t° m > bulk ! the 
reconstruction energy AE^xN (here defined to be posi- 
tive if a reconstruction is favored) is 



5xN 



E^ h N - 5NEf^ -(N + 6)E 



lab 



patom,bulk 



(1) 



For all energies in this work, we chose accurate com- 
putational settings to guarantee a cumulative error of 
Ai?5xAr below 0.02 eV/(lxl) at most. Technical choices 
and convergence tests are summarized in the EPAPS 55 1 
supporting material submitted with this work. 
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FIG. 2: Change of surface energy with reconstruction, com- 
pared to (lxl) [negative of A£ 5x jv, Eq. JTJ], for Pt(100) (up- 
per panel) and Au(100) (lower panel) using DFT-LDA/PBE, 
as a function of the lateral atomic packing density 6. Dash- 
dotted lines: Experimental estimates for the packing den- 
sity from SXRD.[35l. [3(| and for the Pt(100) reconstruction 
energy from the reanalysis of CO titration calorimetry 
data.Hlil 



Figure [5] summarizes our results regarding AEs x n 
for Pt(100) and Au(100), using the approximant cells 
Fig. QJi and b. Since the density of (5xN) planes in- 
creases as TV decreases, we plot AE$ function 
of the lateral atomic packing density relative to the un- 
derlying (100) lattice: 9 — ntex/^ioo- Experimentally, 
SXRD yields 9 wl.242-1.250 for Pt(100)[35|,|3^ and 1.259 
for Au(100).[H, [3^] We find that the two-dimensional 
(5xiV) compression is indeed energetically favored over 
(5x1) both in LDA and PBE. In essence, the energy sur- 
face shows shallow minima that are precisely in line with 
the experimental analysis. The experimental trend of 
a denser "hex" layer for Au than for Pt is clearly con- 
firmed. For Pt, the energy gain of (5xA^) over (5x1) is 
small (~0.01 eV/lxl), but Au roughly doubles its re- 
construction energy. This change brings the reconstruc- 
tion energies of both surfaces close to one another, and 
also to the experimental reconstruction energy estimate 
for Pt(100). In numbers, we find A£:5 X 7v=0.07 eV/lxl 
(Au) vs. 0.10 eV/lxl (Pt) in DFT-LDA. In PBE, the 
reconstructed surface energies are even slightly lower: for 
Pt, AE 5xN =0.11 cV/lxl (PBE) is in remarkably close 
agreement with experiment (0.12±0.02 eV/lxl [la]), 
and we recall that we even expect a further lowering 
by «0.01 eV/lxl for the theoretical value from relax- 
ation beyond the second layer, as indicated by (5x1 ) ap - 



proximants with thicker slabs (see EPAPS material 55]). 
We thus conclude that there is no contradiction between 
theory and experiment for Pt(100), and that the recon- 
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FIG. 3: (a) Top view of the "hex" layer on an Au(100)- 
"(5x20)" reconstruction model (LDA). z coordinates of all 
atoms are color-coded in units of dbuik, relative to the aver- 
age z position of the plane, (b) Height profiles through the 
most corrugated and least corrugated surface areas. 



struction energies of Au(100) and Pt(100) are similar 
when the same in vacuo conditions are applied. In our 
view, the seeming agreement for Au(100)-(5 x 1) theory in 
vacuo and electrochemical experiments was accidental. 

We next consider the energy difference between the 
partially and fully rotated approximants in Figs. [JJd 
and c for Au(100) in LDA at the optimum packing den- 
sity (5x20). The energy gain amounts to 4 meV/lxl 
in favor of the approximant closest to experiment, Fig. 
[TJ;. Between a completely unrotated (5x20) cell and the 
uniformly rotated cell of Fig. [TJ;, the change would be 
8 meV/lxl. This number is still significant when taken 
per reconstructed unit cell, but is also an order of mag- 
nitude smaller than the full reconstruction energy. It is 
thus fully consistent with the observation that the hex 
plane can be modified depending on its environment, 
U H^l without disrupting the reconstruction per se, and 
gives a quantitative idea of the energy scale of these pro- 
cesses. 



In the accompanying EPAPS material, 55] we demon- 
strate the excellent agreement of the overall reconstruc- 
tion geometry calculated here with experimental char- 
acteristics available in the literature. To illustrate the 
consequences for local surface structure variations, Fig. 
[3] shows a top view ("hex" layer only) of the "(5x20)" re- 
constructed Au(100) surface in LDA. Qualitatively, this 
picture looks strikingly similar to published atomically 
resolved STM images. [29| For a more quantitative view, 
the frame below displays the z coordinate of atomic 
zigzag rows in the "5" direction with the largest and 
smallest corrugation amplitude, respectively. The vari- 
ation shows that very different local environments exist 
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Material 


Ir 


Pt 


An 


Au(NR) 


Ag 


Reconstruction 


5x1 


"hex" 


"hex" 


none 


none 


Band center shift [eV] 


-0.19 


-0.24 


-0.13 


-0.01 


-0.08 



TABLE I: Valence band center shift of the projected densities 
of states (LDA) of the surface atoms between (100)-(lxl) and 
(5x1) "hex" approximants, for various elements, compared to 
the computed reconstruction tendency. Au(NR) denotes Au, 
but with a non-relativistic kinetic energy. 



in the surface as a result of the two-dimensional recon- 
struction. The maximum corrugation in the "5" direction 
is 0.65 A (32% of the bulk interlayer spacing dbulk)> com- 
pared to only 15% in the least corrugated area. These 
numbers are very similar in LDA and PBE, and also very 
similar to Pt(100)-(5x25) (29% and 14%, respectively). 
It is precisely such different local environments that drive 
the preferential nucleation and nanostructuring phenom- 
ena on such surfaces. 

Finally, we turn briefly to the electronic-structure 
changes that accompany the reconstruction, which we 
can address directly; such changes are the basis to under- 
stand the role of catalytic activity in 5c? based systems. Q 
For 5d metals, relativity enhances the participation of the 
valence d electrons by shifting them upwards towards the 
s levels. 56, 57] Since d electrons are thus more easily pro- 
moted to sp-like states, enhanced sp bonding has been 
proposed as a cause of the reconstruction. [H, On 
the other hand, enhanced bonding among the d states 



themselves was favored in Ref . [4l| . Table Q] shows the 
shift of the valence band centers of the surface layer pro- 
jected densities of states of Ir(100), Pt(100), and Au(100) 
(all three reconstruct) and Au(100) without relativity, or 
Ag(100) (the latter two surfaces would not reconstruct), 
when going from a (lxl) to a (5x1) reconstruction. In 
short, a characteristic downward shift of the valence band 
center occurs for surfaces that reconstruct, but is much 
smaller for the other two surfaces. An additional angular 
momentum decomposition of the changes shows only a 
tiny transfer of electrons between the d and sp channels, 
so the associated change must happen within the d states 
themselves. Taken together, this analysis thus supports 
the relativistically enhanced d-d hybridization mecha- 
nism suggested by TCH 41 1 over sp promotion. 58, Hj| 



In conclusion, we have demonstrated that large recon- 
struction approximants are indeed necessary to capture 
the subtle reconstruction energy balance of Au(100)- and 
Pt(100)-"hex" quantitatively. Together with the increas- 
ing power of computers and computational methods, re- 
alistically large simulations of chemical processes in the 
presence of reconstruction/deconstruction are now within 
reach. 
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